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Abstract 

Constructive methods for matrices of multihomogeneous (or multigraded) resul- 
tants for unmixed systems have been studied by Weyman, Zelevinsky, Sturmfels, 
Dickenstein and Emiris. We generalize these constructions to mixed systems, whose 
Newton polytopes are scaled copies of one polytope, thus taking a step towards sys- 
tems with arbitrary supports. First, we specify matrices whose determinant equals 
the resultant and characterize the systems that admit such formulae. Bezout-type 
determinantal formulae do not exist, but we describe all possible Sylvester-type and 
hybrid formulae. We establish tight bounds for all corresponding degree vectors, and 
specify domains that will surely contain such vectors; the latter are new even for 
the unmixed case. Second, we make use of multiplication tables and strong duality 
theory to specify resultant matrices explicitly, for a general scaled system, thus in- 
cluding unmixed systems. The encountered matrices are classified; these include a 
new type of Sylvester-type matrix as well as Bezout-type matrices, known as partial 
Bezoutians. Our public-domain Maple implementation includes efficient storage of 
complexes in memory, and construction of resultant matrices. 
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1 Introduction 



Resultants provide efficient ways for studying and solving polynomial systems by 
means of their matrices. They are most efficiently expressed by a generically non- 
singular matrix, whose determinant is a multiple of the resultant, so that the de- 
terminant degree with respect to the coefficients of one polynomial equals that of 
the resultant. For two univariate polynomials there are matrix formulae named af- 
ter Sylvester and Bezout, whose determinant equals the resultant; we refer to them 
as detererminantal formulae. Unfortunately, such determinantal formulae do not 
generally exist for more variables, except for specific cases; this is the topic of our 
paper. 

We consider the sparse (or toric) resultant, which exploits a priori knowledge on the 
support of the equations. Matrix formulae have been studied for systems where the 
variables can be partitioned into groups so that every polynomial is homogeneous 
in each group, i.e. mixed multihomogeneous, or multigraded, systems. This study is 
an intermediate stage from the theory of homogeneous and unmixed multihomo- 
geneous systems, towards fully exploiting arbitrary sparse structure. Multihomoge- 
neous systems are encountered in several areas, e.g. [3,12,8]. Few foundational works 
exist, such as [13], where bigraded systems are analyzed. Our work continues that 
of [7,14,16], where the unmixed case has been treated, and generalizes their results 
to systems whose Newton polytopes are scaled copies of one polytope. These are 
known as generalized unmixed systems, and allow us to take a step towards systems 
with arbitrary supports. This is the first work that treats mixed multihomogeneous 
equations, and provides explicit resultant matrices. 

Sparse resultant matrices are of different types. On the one end of the spectrum 
are the pure Sylvester-type matrices, filled in by polynomial coefficients; such are 
Sylvester's and Macaulay's matrices. On the other end are the pure Bezout-type ma- 
trices, filled in by coefficients of the Bezoutian polynomial. Hybrid matrices contain 
blocks of both pure types. 

We examine Weyman complexes (defined below) , which generalize the Cayley-Koszul 

complex and yield the multihomogeneous resultant as the determinant of a complex. 
These complexes are parameterized by a degree vector m. When the complex has 
two terms, its determinant is that of a matrix expressing the map between these 
terms, and equals the resultant. In this case, there is a determinantal formula, and 
the corresponding vector m is determinantal. The resultant matrix is then said to be 
exact, or optimal, in the sense that there is no extraneous factor in the determinant. 
As is typical in all such approaches, including this paper, the polynomial coefficients 
are assumed to be sufficiently generic for the resultant, as well as any extraneous 
factor, to be nonzero. 

In [16], the unmixed multihomogeneous systems for which a determinantal formula 
exists were classified, but no formula was given; see also [10, Sect. 13. 2]. Identify- 
ing explicitly the corresponding morphisms and the vectors m was the focus of [7]. 
The main result of [14] was to establish that a determinantal formula of Sylvester 
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type exists (for unmixed systems) precisely when the condition of [16] holds on the 
cardinalities of the groups of variables and their degrees. In [14, Thm.2] all such 
formulae arc characterized by showing a bijection with the permutations of the vari- 
able groups and by defining the corresponding vector m. This includes all known 
Sylvester-type formulae, in particular, of linear systems, systems of two univariate 
polynomials, and bihomogeneous systems of 3 polynomials whose resultant is, re- 
spectively, the coefficient determinant, the Sylvester resultant and the classic Dixon 
formula. 

In [14], they characterized all determinantal Cayley-Koszul complexes, which are 
instances of Weyman complexes when all the higher cohomologies vanish. In [7], 
this characterization is extended to the whole class of unmixed Weyman complexes. 
It is also shown that there exists a determinantal pure Bczout-type resultant for- 
mula if and only if there exists such a Sylvester-type formula. Explicit choices of 
determinantal vectors are given for any matrix type, as well as a choice yielding 
pure Bezout type formulae, if one exists. The same work provides tight bounds for 
the coordinates of all possible determinantal vectors and, furthermore, constructs a 
family of (rectangular) pure Sylvester-type formulae among which lies the smallest 
such formula. This paper shall extend these results to unmixed systems with scaled 
supports. 

Studies exist, e.g. [3], for computing hybrid formulae for the resultant in specific 
cases. In [1], the Koszul and Cech cohomologies are studied in the mixed multiho- 
mogeneous case so as to define the resultant in an analogous way to the one used in 
Section 2. In [5], hybrid resultant formulae were proposed in the mixed homogeneous 
case; this work is generalized here to multihomogeneous systems. Similar approaches 
are apphed to Tate complexes [4] to handle mixed systems. 

The main contributions of this paper are as follows: Firstly, we establish the analog 
of the bounds given in [7, Sect. 3]; in so doing, we simplify their proof in the un- 
mixed case. We characterize the scaled systems that admit a determinantal formula, 
either pure or hybrid. If pure determinantal formulae exist, we explicitly provide 
the m- vectors that correspond to them. In the search for determinantal formulae 
we discover box domains that consist of determinantal vectors thus improving the 
wide search for these vectors adopted in [7]. We conjecture that a formula of mini- 
mum dimension can be recovered from the centers of such boxes, analogous to the 
homogeneous case. 

Second, we make the differentials in the Weyman complex explicit and provide de- 
tails of the computation. Note that the actual construction of the matrix, given the 
terms of the complex, is nontrivial. Our study has been motivated by [7], where 
similar ideas are used in the (unmixed) examples of their Section 7, with some con- 
structions which we specify in Example 4.3. Finally, we deliver a complete, publicly 
available Maple package for the computation of multihomogeneous resultant matri- 
ces. Based on the software of [7], it has been enhanced with new functions, including 
some even for the unmixed case, such as the construction of resultant matrices and 
the efficient storage of complexes. 
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The rest of the paper is organized as follows. We start with sparse multihomo- 
geneous resultants and Weyman complexes in Section 2 below. Section 3 presents 
bounds on the coordinates of all determinantal vectors and classifies the systems 
that admit hybrid and pure determinantal formulae; explicit vectors are provided 
for pure formulae and minimum dimension choices are conjectured. In Section 4 we 
construct the actual matrices; we present Sylvester- and Bezout-type constructions 
that also lead to hybrid matrices. We conclude with the presentation of our Maple 
implementation along with examples of its usage. 

Some of these results have appeared in preliminary form in [9] . 



2 Resultants via complexes 

We define the resultant, and connect it to complexes by homological constructions. 
Take the product X := P'^ x ■ ■ ■ x P''' of projective spaces over an algebraically 
closed field F of characteristic zero, for r G N. Its dimension equals the number of 
affine variables n = J2k=i We consider polynomials over X of scaled degree: their 
multidegree is a multiple of a base degree d — (di, . . . , d^) £ N'', say deg/, = Sid. 
We assume sq < ••• < Sn and gcd(so, . . . , s„) = 1, so that the data l,d,s — 
(sq, . . . , Sn) G N"^^ fully characterize the system. We denote by S{d) the vector space 
of multihomogeneous forms of degree d defined over X. These are homogeneous of 
degree dk in the variables Xk ior k = 1, . . . ,r. By a slight abuse of notation, we also 
write S{dk) C P''^ for the subspace of homogeneous polynomials in variables, of 
degree dk- A system of type (Z, d, s) belongs to V — S{sod) ® • • • © S{snd). 

Definition 2.1. Consider a generic scaled multihomogeneous system / = (/o, . . . , /„) 
defined by the cardinahties I e W, base degree d e W and s e N""*"^. The multi- 
homogeneous resultant 7?.(/o, . . . , /„) = T^i,d,sifo, ■ ■ ■ , fn) is the unique up to sign, 
irreducible polynomial of , which vanishes if and only if there exists a common 
root of /o, . . . , /„ in X. 

This polynomial exists for any data I, d, s, since it is an instance of the sparse 
resultant. It is itself multihomogeneous in the coefficients of each /j, with degree 
given by the multihomogeneous Bezout bound: 

Lemma 2.2. The resultant polynomial is homogeneous in the coefficients of each 
fi, i — 0, . . . ,n, with degree 




Proof. The degree degf. TZ of TZ{f) with respect to fi is the coefficient of ■ ■ - ylr 
in the new polynomial: 



JJ(sj(ii7/i H h Sjdryr) = n ^jidiyi H ^ drVr) = -{diyi H h drVrT. 
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In [14, Sect.4] the coefficient of yi ■ ■ ■ yll' in (diyi H h drVrY sfiown to be equal 

to 



ti, . . . , i^y 



tfius proving tfie formula in the unmixed case. Hence the coefficient of y!^ • • • y^^ 
our case is this number multiplied by -^—^ -. □ 



This yields the total degree of the resultant, that is, I]"=o ^eg/- TZ. 

The rest of the section gives details on the underlying theory. The vanishing of the 
multihomogeneous resultant can be expressed as the failure of a complex of sheaves 
to be exact. This allows to construct a class of complexes of finite-dimensional vector 
spaces whose determinant is the resultant polynomial. This definition of the resultant 
was introduced by Cayley [10, App. A], [15]. 

For u & U', H'^ (X, Ox{u)) denotes the q-th cohomology of X with coefficients in 
the sheaf Ox{u). Throughout this paper we write for simplicity H'^{u), even though 
we also keep the reference to the space whenever it is different than X, for example 
i7°(P''=, Mfc). To a polynomial system / = (/o, . . . , fn) over V, we associate a finite 
complex of sheaves K, on X : 

^ ^ • . . ^ ^ ^ • • • ^ ^ (1) 

This complex (whose terms are defined in Definition 2.3 below) is known to be exact 
if and only if /o, . . . , /„ share no zeros in X; it is hence generically exact. When 
passing from the complex of sheaves to a complex of vector spaces there exists a 
degree of freedom, expressed by a vector m = (mi, . . . , rrir) G Z*". For every given / 
we specialize the differentials 6i : Ki ^ -f^j-i, i = l— n, ...,n+lby evaluating at 
/ to get a complex of finite-dimensional vector spaces. The main property is that 
the complex is exact if and only if 7?.(/o, ...,/„) 7^ [15, Prop. 1.2]. 

The main construction that wc study is this complex, which we define in our setting. 
It extends the unmixed case, where for given p the direct sum collapses to (^p) 
copies of a single cohomology group. 

Definition 2.3. For m e 7/', v — —n, . . . ,n-\-l and p = 0, . . . , n -|- 1 set 

0<h<-<ip<n \ e=l ) 

where the direct sum is over all possible indices i\ < ■ ■ ■ < ip. The Weyman complex 

n+l 

K, — K,{1, d, s, m) is generically exact and has terms K^, = K^p. 

p=0 

This generalizes the classic Cayley-Koszul complex. The determinant of the com- 
plex can be expressed as a quotient of products of minors from the 5^. It is invari- 
ant under different choices of m e Z'" and equals the multihomogeneous resultant 

Wo,...,/n). 
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2.1 Combinatorics of K, 



We present a combinatorial description of the terms in our complex, applicable to 
the unmixed case as well. For details on the co-homological tools that we use, see 
[10]. 

By the Kiinneth formula, we have the decomposition 

ifce{ojfc} r 

H^cx)= <S)H^'{^'\ak), (2) 

Hr=q k=l 

where q = p — u and the direct sum runs over all integer sums ji H — ■ + jr = Q, jk G 
{QJk}- In particular, i7°(P'*,Q;fc) is isomorphic to S{ak), the graded piece of P'*^ 
in degree ctfe or, equivalently, the space of all homogeneous polynomials in Zfe + 1 
variables with total degree ctfc, where a — m — zd e 7/" ior z e Z. 

By Serre duality, for any a e Z^, we know that 

Hi{a) ~ - 1 - a)*, (3) 

where * denotes dual, and 1 € N'' a vector full of ones. Therefore H^{ak)* — 
Hh-j^_ak-l-lk). 

Furthermore, we identify (P''^ , a^) as the dual space S{—ak — Ik — I)*- This is the 
space of linear functions A : S{ak) — )■ F. Sometimes we use the negative symmetric 
powers to interpret dual spaces, see also [16, p. 576] . This notion of duality is naturally 
extended to the direct sum of cohomologies: the dual of a direct sum is the direct 
sum of the duals of the summands. The next proposition (Bott's formula) implies 
that this dual space is nontrivial if and only if —ak — ^fe — 1 > 0. 

Proposition 2.4. [2] For any a e Z** and /c G {1, . . . , r}, 

(b) 7^ ^ ait < -k, dimi/''=(r^Q;fc) = (-7^-'). 

(c) //°(P'S ak)y^0^ak> 0, dimi^■°(P'^ ak) = (7^'*)- 

Definition 2.5. Given l,d E W and s E N"+\ define the critical degree vector 
p e N'' by pfe := 4 ELo se - h - 1, for all A; = 1, ... r. 

The Kiinneth formula (2) states that H'^{oi) is a sum of products. We can give a 
better description: 

Lemma 2.6. If H'^[ot) is nonzero, then it is equal to a product H^^ {¥'■'', ak) (8) • • • ® 
i^">(P'^ ak) for some integers ji, . . .,jr with jk E {0, k}, El=iJk = Q- 

Proof. By Proposition 2.4(a), only i/°(P''=, a^) or if''=(P''=, a^) may be nonzero. By 
Proposition 2.4(b,c) at most one of them appears. □ 
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Combining Lemma 2.6 with Definition 2.3 and (2) we get 

<^H^'^(t\mk-jZsi,di\ (4) 

0<ii<---<ip<nfc=l V e=l I 

for some integer sums ii + • • • + = P — 3k £ {0, lk\ such that all the terms in 

r 

the product do not vanish. Consequently, dim (a) = dim {f^^ , a^^ . The 

k=l 

dimension of K^^p follows by taking the sum over all a = m — Y,9=i ^ie^, for all 
combinations {ii < ■■■ <ip} C{0,...,n}. 

Throughout this paper wc denote [u. v] := {u, u + 1, . . . , v}; given p e [0, n + 1], the 
set of possible sums of p coordinates out of vector s is 



=1 



and by convention 5*0 = {0}. By Proposition 2.4, the set of integers z such that both 
i^"°(P'^mfe - zdk) and W'^{¥\mk - zdk) vanish is: 



\ dk^ dk 



nz. 



We adopt notation from [16]: for u e Z, Pk < u «^=^ u > and Pk > u <^==^ 

u < Note that wc use this notation even if Pk = 0. As a result, the z E Z that 
lead to a nonzero H^'' {¥'■'', rrik — zdk), ioi jk = Ik or jk = 0, and p G [0, n + 1], lie in: 

Qp^Sp\ U{Pk, and Q = u;+^Qp. (5) 

Now #Pfc < Ik implies #(UfePfc) < J^- So #{UpSp) > n + 2 implies #(5 > 2. We 
define a function g : Q — > [0, n] by 

?(^) := E (6) 

Pk<z 

Observe that H^{X,m — zd) ^ <^=^ z E Q and j = ^(-z); also the system is 
unmixed if and only if Sp = {p}. Clearly 1 < ^Sp < ("p^), the former inequality 
being strict for s 7^ 1 G N"+^ and p 7^ 0, n + 1. 

The following lemma generalizes [16, Prop. 2. 4]. 

Lemma 2.7. Let G Z, p G {0, . . . , n + 1} and K^^^p given by Definition 2.3; then 
Ku,p 7^ <(=^ G {p - : z G Qp}- 



Proof. Assuming K^^^p ^ 0, there exists a nonzero summand ^(m — zd) 7^ 0. By 
Lemma 2.6 it is equal to i?^i(P''=,mjt-2;4)(8)---(8)i7-''"(P^Smjk-2;4) ^ 0, jk G {0,lk} 
and 

r 

P - = E .^fe = E ^-fc ^ ^ = ~ H ^k- 

k=l Pk<Z Pk<Z 
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Conversely, if z/ G {p — q{z) : z G Qp} then Qp ^ 0. Now z E Qp implies z ^ P, 
which means H'^'^^\m — zd) ^ 0, the latter being a summand of Kj,^p. □ 

One instance of the complexity of the mixed case is that in the unmixed case, given 
p e [0, n + 1], there exists at most one integer v such that K^^p ^ 0. 

All formulae (including determinantal ones) come in dual pairs, thus generalizing 
[7, Prop.4.4]. 

Lemma 2.8. Assume m, m' G Z' satisfy m + m! = p, where p is the critical 
degree vector. Then, K^(m) is dual to Ki_^{'m') for all i/ G Z. In particular, m is 
determinantal if and only if m' is determinantal, yielding matrices of the same size, 
namely dim{Ko{m)) — dim{Ki{m')). 

Proof. Based on the equality m + m' = p we deduce that for all J C [0, n + 1], it 
holds that m' — Y^i^j Sid = —I — 1 — (m — J2iifj Sid). Therefore, for all g = 0, . . . , n, 
Serre's duality (3) implies that H'^{X,m' — J2ieJ ^i^) and m — Sid) 

are dual. 

Let 7^ J = p and u = p — q; since {n + 1 — p) — [n — q) = 1 — {p — q) = 1 — 
we deduce that Ky^p{m) is dual to i^ri_i,^„_|_i_p(m') for all p G [0,^ + 1] which 
leads to Ky{rn)* ~ Xi_^(m') for all G Z, as desired. In particular, K_i{ra) ~ 
K2{m!) and KQ{m) ~ Kl{m'), the latter giving the matrix dimension in the case 
of determinantal formulae. □ 



3 Determinantal formulae 

This section focuses on formulae that yield square matrices expressing the resul- 
tant without extraneous factors and prescribes the corresponding determinantal m- 
vectors. 

Determinantal formulae occur only if there is exactly one nonzero differential, so 
the complex consists of two consecutive nonzero terms. The determinant of the 
complex is the determinant of this differential. We now specify this differential; for 
the unmixed case see [16, Lem.3.3]. 

Lemma 3.1. If m G Z'' is determinantal then the nonzero part of the complex is 

Proof. The condition that m is determinantal is equivalent to the fact that N :— 
{p — q{z) : z & Qp,p & [0,n + 1]} consists of two consecutive integers. 

Let Zi = minQ < Z2 — maxQ, since > 2. There exist pi,P2 with pi < p2 such 
that Zi — Z^^Li and Z2 — Sa=i where the indices are sub-sequences of [0, n], 



8 



of length pi and p2 resp. The pi integers 

0, -So, So + Si, . . . , So + • • • + •Spj_2 G Z 

arc distinct, smaller than zi, hence belong to Upi^^ziPk- Also, it is clear that, for all 
ke[l,r], #Pfc < f/fe/4] < Ik thus 

Pi < # U Pk< Y.*Pk< q{zi). (7) 

Pk<zi Pk<zi 

This means pi — q{zi) < 0. Similarly, the n + 1 — p2 integers s„ H — • + sq, . . . , s„ + 
• • • + s„_p2 are distinct, larger than Z2, hence belong to Up^^zi^k, so: 

n-p2 + l<# U Pk< *Pk< Y.^k- (8) 

Pk>Z2 Pk>Z2 Pk>Z2 

This means n + 1— p2 <n — q{z2), thus p2 — 5(^2) > 1- Hence there exists a positive 
integer in N; from (7) we must have a non-positive integer in N. Since 771^ A?^ = 2 and 
the integers of N are consecutive we deduce that N — {0,1}. □ 

Corollary 3.2. If m G Z*" is determinantal, then equality holds in (7). In particular, 
for k e [1, r] s.t. either < Z\ or > Z2, we have #Pfe = and any two such P^ 
are disjoint. 

Proof. Lemma 3.1 combined with Lemma 2.7 imply Pi — q{zi) > 0, and (7) im- 
plies pi — q{zi) < 0, hence we deduce pi — q{zi) = 0. Now equality in (7) gives 
J2p^<zi #Pk = q{zi) = Ep,<2i h; combining with #Pfc < k we deduce #Pfe = k for 
all k in this sum. Similarly for n -|- 1 — P2 = 'E,p^.>z2 h = ^ ~ ?(^2)- CH 

3. 1 Bounds for determinantal vectors 

We generahze the bounds in [7, Sect. 3] to the mixed case, for the coordinates of all 
determinantal m— vectors. We follow a simpler and more direct approach based on 
a global view of determinantal complexes. 

Lemma 3.3. If a vector m e Z'' is determinantal then the corresponding (Ji Pk is 
contained in [0,X^QSj]. 

Proof. It is enough to establish that P^ > and P^ < Z^o for all /c G [1,'r]. Proof 
by contradiction: Let m be a determinantal vector, and P = UgPfc. Let Zi.,Z2 as in 
the proof of Lemma 3.1. If Zi < P^ < Z2 it is clear that < Zi < Pk < Z2 < Z^q ^ 

If Z2 < Pk, Corollary 3.2 implies #Upj.>22-Pfe — H^P^ where R := {s„ + ■ ■ • -|- 
So; • • • 5 "S„ -|- ■ ■ ■ -|- Sfj-pj}. By the definition of Z2, R C Upj.>z2 Pk, thus Upj.>z2 Pk = 

R ^ [0,Eo Si]- Similarly, Up^<2i Pk = {O,so,so + Si, . . . ,soH hSpi-2} ^ [0,Eo s^], 

which proves the lemma for Pk < zi. □ 
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The bound below is proved in [7, Cor. 3. 9] for the unmixed case. They also show 
with an example that this bound is tight with respect to individual coordinates. We 
give an independent, significantly simplified proof, which extends that result to the 
scaled case. 

Theorem 3.4. For determinantal m e Z'', for all k we have 

n 

max{-dk, -Ik} <mk<dk^Si-l + min{4 - Ik, 0}. 





Proof. Observe that by Lemma 3.1 there are no A; e [l,r] such that < or 
Pk > J2q Si- Combining this fact with Lemma 3.3, we get 

n 

rrik/dk > -1 and {nik + lk)/dk <l + ^Si (9) 



for all k e [1, r]. Furthermore, the sets P^, A; e [1, r] can be partitioned into two (not 
necessarily non-empty) classes, by considering the integers Zi, Z2 of Lemma 3.1: 



• Pk < Zi or Pk > Z2, with cardinalities #Pfc = Ik- 

• zi < Pk < Z2, without cardinality restrictions (possibly empty). 



Taking into account that Pj^ = ^^^j n Z we get 



{ruk + lk)/dk > and mfc/4 <^Si 



(10) 



for all k e [1, r]. 



□ 



Our implementation in Section 5 conducts a search in the box defined by the above 
bounds. For each m in the box, the dimension of K2 and K^i is calculated; if both 
are zero the vector is determinantal. Finding these dimensions is time consuming; 
the following lemma provides a cheap necessary condition to check before calculating 
them. 

Lemma 3.5. If m G Z^ is determinantal then there exist indices k,k' e [l,r] such 
that rUk < dk{sn-i + Sn) and rUk' > dk' Eo~^ - k'- 



Proof. If for all /c, ruk/dk > Sq + Si then q{sn-i + s„) = by (5), so for p = 2 we 
have p — g'(s„_i -|-s„) = 2 — = 2 which contradicts the fact that m is determinantal. 
Similarly, if for all k, (ruk + lk)/dk < So~^ g — n and for p = n — 1 

we have p — q {^2 Si) — {n — 1) — n — —1, which is again infeasible. □ 
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3.2 Characterization and explicit vectors 

A formula is determinantal if and only if K2 — K_i = 0. In this section we provide 
necessary and sufficient conditions for the data l,d,s to admit a determinantal 
formula; we call this data determinantal. Also, we derive multidimensional integer 
intervals (boxes) that yield determinantal formulae and conjecture that minimum 
dimension formulae appear near the center of these intervals. 

Lemma 3.6. If m e Z'' is a determinantal vector for the data I, d, s, then this data 
admits a determinantal vector m' e Z'' with f] Pk' — for all k,k' e [1, r]. 

Proof. Suppose nii/di < ruj/dj. Let Pi{m) nPj{m) = [u,v] C Z. Set m'j — mj +tdj 
where i e Z is the minimum shift so that Pi{m') fl Pj{m') — and Pj{m') satisfies 
Theorem 3.4. For all k ^ j, let m'^ — m^. 

Any vector in 7/ defines a nontrivial complex, since Q ^ 0. \a particular, m' 
is determinantal because P{m) C P{m!), i.e. no new terms are introduced, but 
possibly some terms vanish. Repeat until all fl Pfe/ — 0. □ 

Let a : [l,r] [l,r] be any permutation. One can identify at most r! classes of 
determinantal complexes, indexed by the permutations of {1, . . . , r}. This classifica- 
tion arises if we look at the nonzero terms that can occur in the complex, provided 
that the sets Pk satisfy 

P.(l) < Pa{2) <■■■< Pair) 

where we set Pj < Pj <^=^ iTii/di < rrij/dj. Any given m defines these sets, 
as well as an ordering between them. This fact allows us to classify determinantal 
m-vectors and the underlying complexes. 

For this configuration, expressed by a, the only nonzero summands of can be 

Ku,v+q where q takes values in the set {0, la{i), la{i) + la{2), .... n}. To see this, observe 
that q{z) = J2p^<z h, z cannot attain more than r + 1 distinct values; so if the 
relative ordering of the P^ is fixed as above, then these are the only possible values 
of q. This leads us to the following description of K2 and K_i: 

k=l k=l 

As a side remark, note that the proof of Lemma 2.8 implies that the dual of K^[m) 
is Kl_^{p — m) where r is the permutation s.t. r(z) := r + 1 — ^{i). 

Let n[k] := J2TT{i)<TT{k) h- If tt = Id this is ld[k] — li + • • ■ + l^- We now characterize 
determinantal data: 

Theorem 3.7. The data Z, cZ, s admit a determinantal formula if and only if there 
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exists TT : [l,r] — >■ [l,r] s.t. 

n 7r[fe-l]+l 

dk XI Si-lk< dk ^ Si, V/c. 

n— 7r[fc]+2 

Proof. We assume without loss of generality that tt = Id. This is not restrictive, since 
if TT 7^ Id then we can re- number the variable groups such that k' :— Tr~^{k). Hence 

if we set Lk := ifo^'"'^^^^ Si = min 5'id[fe-i]+2 and Rk := En-id[fe]+2 = max5'id[fe]-i 
then the relations become: 

dkRk — h < dkLk, Vfc. 

Throughout this proof, whenever we use non-positive indices j < for Ij or Id[j], 
these quantities will be zero, and the results in this case are straightforward to 
verify. Also, note that the dual complex is given by the "reversed" permutation, and 
in particular, Kl^ ~ K2, therefore any results on the nullity of K^i can be directly 
used to prove the nullity of K2. 

(^) Assume that the inequalities hold. Then for all k there exists an integer rrik 
such that 

dkRk - lk<fnk< dkLk - 1. (12) 
Let m = (mi, . . . ,771^). We shall prove that this vector gives a determinantal for- 
mula; it suffices to show that for all k e [l,r], i^2,2+id[fe-i] = -f^-i,-i+id[ifc] = 0, since 
in (11) we have a = 7r~^ = Id. 

• If Zfe > 3, we have 

Id[A;] - 1 = Id[A; - 1] + - 1 > Id[A; - 1] + 2. (13) 



Thus Lk < Rk, since Sj > 1; also our hypothesis (12) translates into the inclusion 

'mk nik + lk 



Lk, Rk 



C 



Pk. 



(14) 



djfc ' dk 

Now, by (13) we derive 

min5'id[fe_i]+2 ^ Lk < min5'id[fc]-i and maxS'id[fe]_i ^ Rk > niax5'id[fe-i]+2 

so (14) imphes S'id[fe-i]+2 ^ as weU as Su[k]-i Q Pk and thus K2,2+id[k-i] = 0, 
X_i _i+id[fe] = by Proposition 2.4. 

• lilk < 2, it is Id[A:] - 1 = Id[A: - 1] + 4 - 1 < ld[k - I] + 1. In this case we wiU 
prove i^2,2+id[fc-i] = 0, _i+id[fc] = using Lemma 2.7. 

Let z E Qp for p = —1 + ld[k]. Prom Rk < (m^ -|- lk)/dk it is clear that Pk ^ thus 
q{z) <Id[A;-l]. Also, 



Id[A; - 2] + 2 = Id[A;] - h-i -lk + '2< ld[k] - 1. 
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where the last inequahty is taken under the assumption max{Zfe, Ik-i} > 2. We treat 
the case — Ik-i — 1 seperately. Hence 

z > min 5'_i+id[fc] > minS'2+id[fc-2] = -^fc-i > ruk-i/dk-i. 

This imphes Pk-i < z ^ q{z) > ld[k — 1]. We conclude that q{z) — ld[k — 1] and 
thus 

p - q{z) = (-1 + Id[A;]) - Id[A; - 1] = -1 + ^ G [0, 1]. (15) 
By Lemma 2.7 we see that _i_|_id[fc] = 0, since p — q{z) ^ —1. 

To complete the proof, suppose Ik — Ik-i — 1 and p = ld[k] — 1. We get Id [A;] — 1 = 
Id[A; — 2] + 1 and therefore z > Pk-2 =^ > Id[A; — 2]. Recall that q{z) is also 
upper bounded by ld[A; — 1]. We derive that for z E Qp it holds p — 1 < q{z) < p, 
therefore p — q{z) G [0, 1], and again by Lemma 2.7, i^_i,_i_)_id[fc] = 0. 

As already pointed out, by using duality one can see that, for z' G Q2+id[k-i], it 
holds p - q{z') ^ 2, therefore -f^2,2+id[fe-i] = 0. 

(=^) Suppose that m G Z'' is determinantal, namely K2{ra) — K_i{m) — 0. 
Lemma 3.6 implies that we may assume the sets Pj are pairwise disjoint. By a 
permutation of the variable groups we also assume that the Pj sets induced by m 
satisfy 

Pl<P2<---<Pr- 

The sets Pj have to be distributed along / := [0, X^o (Lemma 3.3) and the con- 
nected components of I\^Pj are subsets of SpUSpj^i, p G [0, n + 1] since they define 
a determinantal complex. In particular, for p — ld.[k] — 1 we get Pk-i < SpU Sp+i < 
Pk+i- Now the definition of Rk as an element of <S'id[fe]-i implies Pk-i < Rk < Pk+i, 
i.e. we have the implications (similarly for L^): 

RkeU[Pj ^ Rk&Pk and LkeU\Pj ^ Lk e Pk- (16) 



Suppose nik < dkRk — h, or equivalently ; < Rk- Then Rk ft Pk, hence by 

dk 

(16) we must have Rk ^ UjPj which leads to z — Rk E Qp, p — ld[k] — 1. This 
implies 

q{z) > ld[k] =^ p- q{z) < ld[k] - 1 - Id[k] = -1 ^ K_i ^ 0, 

which is a contradiction. In the same spirit, if > dkLk, we are led to z' — Lk & Qp, 
p = Id[A; - 1] + 2, then 

q(z') < Id[k - 1] => p- q{z') >p- Id[k - 1] = 2 K2 0, 

which again contradicts our hypothesis on m. 

We conclude that any coordinate of m must satisfy dkRk ~ h ^ < dkLk-, 
hence the existence of m implies the inequality relations we had to prove. □ 
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Corollary 3.8. For any permutation tt : [l,r] — >■ [l,r], the vectors m e con- 
tained in the box 



n 



7r[A;-l]+l 



dk J2 Si-lk<mk<dk Yl 

n— 7r[fe]+2 



for k = 1, . . . , r are determinantal. 

It would be good to have a characterization that does not depend on the per- 
mutations of [l,r]; this would further reduce the time needed to check if some 
given data is determinantal. One can see that if r < 2 an equivalent condition 

n 

is dk ^ Si — Ik < dk{so + Si) for all k e see [5, Lem.5.3] for the case 

r = 1. It turns out that for any r E N this condition is necessary for the existence 
of determinantal vectors, but not always sufficient: the smallest counterexample is 
I — (1,2,2), d = (1,1,1), s = (1,1,1,1,2,3): this data is not determinantal, al- 
though the condition holds. In our implementation this condition is used as a filter 
when checking if some data is determinantal. Also, [5, Cor. 5. 5] applies coordinate- 
wise: if for some k, Ik > 7 then a determinantal formula cannot possibly exist unless 
dk — 1 and all the Sj's equal 1, or at most, = s„ = 2, or all of them equal 1 
except s„ = 3. 

We deduce that there exist at most r! boxes, defined by the above inequalities that 
consist of determinantal vectors, or at most r!/2 matrices up to transpose. One can 
find examples of data with any even number of nonempty boxes, but by Theorem 3.7 
there exists at least one that is nonempty. 

If r = 1 then a minimum dimension formula lies in the center of an interval [5] . We 

conjecture that a similar explicit choice also exists for r > 1. Experimental results 
indicate that minimum dimension formulae tend to appear near the center of the 
nonempty boxes: 

Conjecture 3.9. If the data l,d,s is determinantal then determinantal degree 
vectors of minimum matrix dimension lie close to the center of the nonempty boxes 
of Corollary 3.8. 

We conclude this section by treating the homogeneous case, as an example. 

Example 3.10. The case r — 1, arbitrary degree, has been studied in [5]. We shall 
formulate the problem in our setting and provide independent proofs. Let n, G Z, 
s e Z"o^. This data define a scaled homogeneous system in P"; given m G Z, we 



obtain P 



m m+n 



n Z. In this case there exist only zero and nth cohomologies; 



zero cohomologies can exist only for i/ > and nth cohomologies can exist only for 
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u < 1. Thus in principle both of them exist for i/ e {0, 1}. Hence, 

K^^^, l<u<n+l 

Ki,,n+i^, -n < i/ < 

i.e. the complex is of the form: 

Kn+l,n+l y Ki^l © i^l,n+l ^ ^^0,0 © i^O,n i^-l,n-l > i^-n,0 ^ 

We can explicitly give all determinantal integers in this case: 

^2 = <^ ^2,2 = <^ Q2^0 <^ S2CP, 

thus 

in in 
min5'2 > — <^=^ sq + si > — <^==^ m < (sq + si)d. 
a a 

Similarly K^i — <^=^ Qn-i — ^ Sn~i Q P and thus 

m + n " 
max < — - — <^=^ m > a} ^Sj — n. 

1=2 

Consequently, a determinantal formula exists iff dYj2 Si — n< (so + si)<^, also verified 
by Theorem 3.7. In this case the integers contained in the interval 



d^Si-n-l , d{sQ + si) j 

1=2 ^ 



are the only determinantal vectors, also verifying Corollary 3.8. Notice that the sum 
of the two endpoints is exactly the critical degree p. 

In [5, Cor. 4. 2, Prop. 5. 6] it is proved that the minimum- dimension determinantal for- 
mula is attained at m = Lp/2J and m — [p/2] , ie. the center(s) of this interval. For 
an illustration see Ex. 4.5. □ 



3.3 Pure formulae 



A determinantal formula is pure if it is of the form Ki^a — ^ -ft'o,b for a, 6 e [0, n + 1] 
with a > b. These formulae are either Sylvester- or Bezout-type, named after the 
matrices for the resultant of two univariate polynomials. 

In the unmixed case both kinds of pure formulae exist exactly when for all G [1, r] 
it holds that min{/fc, d^} = 1 [14,7]. The following theorem extends this characteriza- 
tion to the scaled case, by showing that only pure Sylvester formulae are possible and 
the only data that admit such formulae are univariate and bivariate-bihomogeneous 
systems. 
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Theorem 3.11. If s ^ 1 a pure Sylvester formula exists if and only if r < 2 and 
Z = (1) or Z = (1, 1). If Zi = n = 1 the degree vectors are given by 

1 

m — di^^Si — 1 and m' — —1, 



whereas if Z = (1, 1) the vectors are given by 

Pure Bezout determinantal formulae cannot exist. 
Notice the duality m -\- m! — p. 



Proof. It is enough to see that if a pure formula is determinantal the following 
inequalities hold 

which implies that equalities hold. The inequality on the left follows from the fact 
that every S'p, p e [0,n + 1] contains at least one distinct integer since the se- 
quence 0, So, So + si, ■ ■ ■ , Ho Si is strictly increasing. For the right inequality, note 
that the vanishing of all K,,^p with p ^ a,b implies Qp = (see Lemma 2.7). Thus 
^p^a,bSp C Ul^iPk so the cardinality is bounded by ^U[Pk < #Pfc < Yli h = n. 
Consequently # ^p^a,b Sp = n. Suppose n > 2; the fact #(S'j U Sj) > 2 for 
{i,j} {0, 1} implies Up^L^^bSp — SiU Sj for some i.e. #{a, 6} = n, contra- 
diction. Thus n <2. 

Take n = 2. Since #(5*0 U 5*3) = 2, the above condition is satisfied for a = 2, 6 = 1: 
it is enough to set U^Pk — SqU S3 = {0,J2o Si}, thus the integers of U^Pfc are not 
consecutive, so r > 1 and I — (1, 1). Similarly, if n = Z = 1 two formulae are possible; 
for U{Pk = -5o = {0} (a = 2, & = 1) or U{Pk ^82^ {sq + (a = 1, 6 = 0). 

All stated m- vectors follow easily in both cases from (m^ + lk)/dk = and {rrik + 
h)/dk — J2o Si. A pure Bezout determinantal formula comes from iCi^n+i — )■ ^^0,0- 
Now UfePfe contains Si \J ■■■ Li Sn hence # Ufe > n. Thus it cannot exist for 

All pure formulae above are of Sylvester-type, made explicit in Section 4. If n = 1, 
both formulae correspond to the classical Sylvester matrix. 

If s = 1 pure determinantal formulae are possible for arbitrary n, r and a pure 
formula exists if and only if for all k, Ik = 1 or dk = 1 [7, Thm.4.5]; if a pure 
Sylvester formula exists for a,b = a — 1 then another exists for a = 1,6 = [7, 
p. 15]. Observe in the proof above that this is not the case if s 7^ 1, n = 2, thus the 
construction of the corresponding matrices for a 7^ 1 now becomes important and 
highly nontrivial, in contrast to [7]. 
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4 Explicit matrix construction 

In this section we provide algorithms for the construction of the resultant matrix 
expressed as the matrix of the differential 5i in the natural monomial basis and we 
clarify all the different morphisms that may be encountered. 

Before we continue, let us justify the necessity of our matrices, using I — d — {1,1) 
and s — (1, 1,2), that is, the system of two bilinear and one biquadratic equation 
to be examined in Example 5.1. It turns out that a (hybrid) resultant matrix of 
minimum dimension is of size 4x4. The standard Bezout-Dixon construction has 
size 6x6 but its determinant is identically zero, hence it does not express the 
resultant of the system. 

The matrices constructed are unique up to row and column operations, reflecting the 
fact that monomial bases may be considered with a variety of different orderings. 
The cases of pure Sylvester or pure Bezout matrix can be seen as a special case 
of the (generally hybrid, consisting of several blocks) matrix we construct in this 
section. 

In order to construct a resultant matrix we must find the matrix of the linear map 
Si : Ki ^ Kq in some basis, typically the natural monomial basis, provided that 
K_i = 0. In this case we have a generically surjective map with a maximal minor 
divisible by the sparse resultant. If additionally K2 = then dimT^i = dimi^o and 
the determinant of the square matrix is equal to the resultant, i.e. the formula is 
determinantal. We consider restrictions : i^i,a ~^ -^0,6 ioi any direct summand 
Ki^a, -K^o,6 of Ki, Ko respectively. Every such restriction yields a block of the final 
matrix of size defined by the corresponding dimensions. Throughout this section the 
symbols a and b will refer to these indices. 

4-1 Sylvester blocks 

The Sylvester- type formulae we consider generahze the classical univariate Sylvester 
matrix and the multigraded Sylvester matrices of [14] by introducing multiplication 
matrices with block structure. Even though these Koszul morphisms are known to 
correspond to some Sylvester blocks since [16] (see Proposition 4.1 below), the exact 
interpretation of the morphisms into matrix formulae had not been made explicit 
until now. We also rectify the Sylvester-type matrix presented in [7, Sect. 7.1]. 

By [16, Prop. 2. 5, Prop. 2. 6] we have the following 

Proposition 4.1. [16] If a — 1 < 6 then 6a,b = 0. Moreover, if a — 1 = 6 then 6a,b is 
a Sylvester map. 

If a = 1 and 6 = then every coordinate of m is non-negative and there are only zero 
cohomologies involved in K^^i = 0j H^{m — Sid) and ^^0,0 = H^{m). This map is a 
well known Sylvester map expressing the multiplication [qq, . . . ,gn) 1 — > Y^Z=o9ifi- 
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The entries of the matrix are indexed by the exponents of the basis monomials of 
0j S{m — Sid) and S{m) as well as the chosen polynomial fi. The entry indexed 
(i, a), /3 can be computed as: 

coef (/j, cc^"") , i = 0, 

where a;" and run through the corresponding monomial bases. The entry (?', a), (3 
is zero if the support of fk does not contain (3 — ol. Also, by Serre duality a block 
Ki^n+i ^ -^^o,n corresponds to the dual of Ki^i -^o,0; i-e. to the degree vector 
p — m, and yields the same matrix transposed. 

The following theorem constructs corresponding Sylvester-type matrix in the general 
case. 

Theorem 4.2. The entry of the transposed matrix of 6a,b '■ -f^i,a -f^o,a-i in row 
(/, a) and column (J, y9) is 

jo, if J ^ /, 

|(_l)fc+icoef(/,„a:-), if /\ J=M, 

where 7 = {ii < i2 < • • • < ia} and J = {ji < < • • • < ja-i}, I, J C {0, . . . , n}. 
Moreover, a,/3 G run through the exponents of monomial bases of H"'^^{m — 
dELi H^-\m - dYXX Sje), and u e N", with Ut = |A - ^tl- 

Proof. Consider a basis of A"^! {^ii,i2....,ia '■ < ii < 22 < ■ ■ ■ < ?'a and 
similarly for tC'^V , where Cq, . . . , e„ is a basis for V . This differential expresses a 
classic Koszul map 

n 

dai^in ■ ■ ■ ■) ^ia) — J^l"!) fik^ii,...,ik-i,ik+i,—,ia 

and by [16, Prop. 2. 6], this is identified as multiplication by /j^, when passing to the 
complex of modules. 

Now fix two sets / C J with I\J = {ik}, corresponding to a choice of basis elements 
e/, cj of the exterior algebra; then the part of the Koszul map from e/ to cj gives 

(-l)'=+iM(/,J : H''-\m - d^s,) ^ H''-\m - dj^se) 

This multiplication map is a product of homogeneous multiplication operators in the 
symmetric power basis. This includes operators between negative symmetric powers, 
where multiplication is expressed by applying the element of the dual space to /j^ . 

To see this, consider basis elements w°', that index a row and column resp. 
of the matrix of M{fi^). Here the part Wk of w associated with the fc-th variable 
group is either x'^*' or a dual element indexed by a*;. We identify dual elements 
with the negative symmetric powers, thus this can be thought as x'^"''' . This defines 
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a./(3 & Z"; the generalized multihomogencous multiplication by fi,^ as in [16, p. 577] 
is, in terms of multidegrees, incrementing \dcj:\ by Sid^ to obtain 1/3^1, and hence 
the corresponding matrix has entry coef (/j^,, a;'"), where Ut = \(3t — (^t\, t G [l,n]. 
The absolute value is needed because for multiplication in dual spaces, the degrees 
satisfy -\ak\ + Sidk = =^ Sidk = la^l - |/3fc| = -{\/3k\ - □ 



In [7, Sect. 7.1], an example is studied that admits a Sylvester formula with a — 
2,b = 1. The matrix derived by such a complex is described by Theorem 4.2 above 
and does not coincide with the matrix given there. The following example is taken 
from there and presents the correct formula. 

Example 4.3. Consider the unmixed case I = (1, 1), d = (1, 1), as in [7, Sect. 7.1]. 
This is a system of three bi-linear forms in two affine variables. The vector m = 

(2,-1) gives i^i = Ki^2 = H\0,-3fi) and Ko = Ko,i = H\l,-2)(^^. The 
Sylvester map represented here is 

^1 ■ {go, 9i, 92) H> (-50/1 - ^1/2, 5o/o - 52/2, ^1/0 + 52/1) 

and is similar to the one in [6]. By Theorem 4.2, it yields the following (transposed) 
matrix, given in 2 x 2 block format: 

-M(/i) M(/o) 
-M{h) M(/o) 
-M{h)M{h) 

If = Co + CiXi + C2X2 + C3X1X2 the matrix of the multiplication map 
M{g) : S{Q)®S{iy3w ^ wg&S{l)®S{Qy 



in the natural monomial basis is 



C2 C3 



Co Ci 



as one can easily verify by hand calculations 



or using procedure multmap of our MAPLE package presented in Section 5. 



□ 



4.2 Bezout blocks 



A Bezout-type block comes from a map of the form 5a,b '■ -f^i,a — > -f^o,6 with a — 1 > 6. 
In the case a — n+1, b = this is a map corresponding to the Bezoutian of 
the system, whereas in other cases some Bezout-like matrices occur, from square 
subsystems obtained by hiding certain variables. 

Consider the Bezoutian, or Morley form [11], of /q, . . . , /„. This is a polynomial of 
multidegree (p, p) in F[i, y] and can be decomposed as 

pi Pr 

A:= ^•••^A„(:r)-r 

U1=0 Ur=0 
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where Au{x) e -S" has degA„(x) — p — u. Here x = {xi, . . . ,Xr) is the set of 
homogeneous variable groups and y — {yi, . . . , y^) a set of new variables with the 
same cardinalities. 

The Bezoutian gives a linear map 

A"+V^ S{p-m)®S{m). 

mk<Pk 

where the space on the left is the (n + l)-th exterior algebra of ^ = S{sQd) © 
• • • ® S{snd) and the direct sum runs over all vectors m e Z'' with < pk for all 
k e [1, r]. In particular, the graded piece of A in degree (p — m, m) in (x, y) is 



^p—Tn,m ■ — ^ ^ Au(!k) ■ y 

Uk=mk 



for all monomials y^ of degree m and coefficients in ¥[x] of degree /o — m. It yields 
a map 

S{p - m)* S{m) 

known as the Bezoutian in degree m of /o, . . . , /„. The differential of Ki^^+i -f^o.o 
can be chosen to be exactly this map, since evidently ^0,0 = H^{m) ~ S{m) and 

according to Serre duality (sec Section 2.1). Thus, substituting the critical degree 
vector, we get = S{p — m)* . 

The polynomial A defined above has n + r homogeneous variables and its homoge- 
neous parts can be computed using a determinant construction in [1], which we adopt 
here. We recursively consider, for A; = 1, . . . , r the uniquely defined polynomials //J^ 
where < j < as follows: 

fi = xi^ofiS + ■ ■ ■ + xi,hfi% fid e F [Xij, Xi,zJ [x2,...,Xr], (17) 

for all z = 1, . . . , n. To define f-'^\ for 2 < k < r and < j < Ij, we decompose 
fih-i '^ith respect to the group Xf 

fij ^ ^ • • • ' ^k-l,lk-i] [Xk,j, Xk,lk] [Xk+1, ...,Xr]. 

Overall we obtain a decomposition 

fi = E xijfS + xi,i, E /g^ + • • • + n E ^k,jfS^ + ■■■ 

j=0 j=0 t=l j=l 

r—l lr — 1 r 

• • • + n ^t,it E ^r,jfij + n ^t,hfiifi,ir e F[xi,i,, ■ ■ ■ , 2:1,;,] 
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of the polynomial fi, for all i = 1, . . . , n. The order of the variable groups, from left 
to right, corresponds to choosing the permutation tt = Id. The determinant of size 
(n + 1) X (n + 1) given by 



V 



Jo,o ■ ■ ■ Jo,h-i ■ ■ ■ Jo,o ■ ■ ■ Jo,ik-i ■ ■ ■ -lOfi ■ ■ ■ JOy 



(k) 



Ar) 



p{l) p(L) r\Kj r\K) pyr) p 

Ji,0 ■ ■ ■ Ji,li-1 ■ ■ ■ Ji,0 ■ ■ ■ Ji,lk-i ■ ■ ■ Jifi ■ ■ ■ J 



dk) 



(fc) 



Ar) 



Ar) 
<i,lr 



M) Ai) Ak) 
J n,0 ■ • • J n,li — l ■ ■ ■ J n,0 



r{k) pyij ryi) 

Jn,/fc — 1 ■ ■ ■ J n,0 ■ ■ ■ J n,l 



(r) 



Ar) 



is equal to Ap_^ „i, in our setting, as we have the following: 

Theorem 4.4. [1] The determinant T> is an inertia form of degree pk — ruk with 
respect to the variable group Xk, k — 1, . . . ,r. 

Let us show a more simple construction of some part Ap-rn,m using an afiine 
Bezoutian. Let Xk = {xk,i, ■ ■ ■ ,Xk,ik) the (dehomogenized) k-th variable group, and 
Uk — {Vk,!, ■ ■ ■ 1 yk,ik)- As ^ result the totality of variables is a; = {xi, . . . , Xr) and 
2/ = (2/l,■■■,yr)■ 
We set Wt, t = 1, . . .n — 1 the conjunction of the first t variables of y and the last 
n — t variables of x. 

lia — n+l,b — the affine Bezoutian construction follows from the expansion of 



foix) foiwi) ■■■ fo{wn-i) My) 



r Ik 

I n \[i.^kj-ykj) 

k=lj=l 



as a polynomial in F[y] with coefficients in F[a;]. Hence the entry indexed q;,/3 of 
the Bezoutian in some degree can be computed as the coefficient of x'^y^ of this 
polynomial. 

We propose generalizations of this construction for arbitrary a, h that are called 
partial Bezoutians, as in [7]. It is clear that a — 1 = q{zi) and b = q{z2), for Zi G Qa 
and Z2 G Qb- The difference a — h — 1 = 1^0=1 he where fci, . . . , /cj is a subsequence 
of [1, r], since if < 6 then Pk < a thus 

q{a) - q{b) = J2 k - J2 ^ = J2 h- 



Pk<a 



Pk<b 



b<Pk<a 



These indices suggest the variable groups that we should substitute in the partial 
Bezoutian. Note that in the case of Bezoutian blocks, it holds a — 6 — 1 > thus 
some substitutions will actually take place. Let ii,...,ia-b be a subsequence of 
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[0,n]. We can define a partial Bezoutian polynomial with respect to fi^, . . ■ , fi^_ 
and yf,^,...,yk^ as 



(18) 



=ii=i 



In this Bezoutian, only the indicated t/-variable substitutions take place, in succes- 
sive columns: The variable vector w differs from x at a^fci, • • • , a^fct, these have been 
substituted gradually with y^^ , • • • , 2/^^ • Note that w generalizes the vectors Wt de- 
fined earlier, in the sense that the variables of only specific groups are substituted. 
The total number of substituted variables is a — 6 — 1, so this is indeed a Bezout 
type determinant. 

For given a and 6, there exist (^i*!^) Partial Bezoutian polynomials. The columns of 
the final matrix are indexed by the cc-part of their support, and the rows are indexed 
by the y-part as well as the chosen polynomials /j^, . . . , 

Example 4.5. Consider the unmixed data 1 = 2, d = 2, s = (1, 1, 1). Determinantal 
formulae are m e [0,3], which is just m = 0, m = 1 and their transposes. Notice 
how these formulae correspond to the decompositions ofp = 3 = 3-|-0 = 2-|-l. 
In both cases the complex is of block type Ki^^ — )■ Kq^ © Ko^2- The Sylvester part 
Ki 3 Kq^2 can be retrieved as in Ex. 4.3. For m = the Bezout part is H^{—6) ~ 
S'(3)* H^{0) ~ 'S'(O), whose 5x1 matrix is in terms of brackets 



[142] [234] + [152] [235] [042] [052] 



A bracket [ijk] is defined as 



[ijk] := det 



Oj aj Ofc 
bi bj bk 

Cj Cj Cj^ 



where ai,bi,Ci denote coefficients of /o,/i,/2 respectively, for instance /2 = cq + 
C1X2 + 02X2^ + C3X1 + + c^xx' . Now, for m = 1 we have Bezout part 5) ~ 

,5(2)* H°{1) ~ S{1), which yields the 5 x 3 matrix 



[142] [152] + [234] [235] [042] [052] 

[152] [154] + [235] [354] [052] [054] 
[132] + [042] [052] + [134] [135] [041] + [032] [051] 



T T 



□ 
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rontinp 


fiinrt.ion 


Make system 


output polynomials of type (/, s) 


mBezou't 


compute the m-Bezout bound 


cL±±jjex V ecs 


enumerate cill determincintcLl tn — vectors 


detboxes 


output tlie vector boxes of Cor. 3.8 


-F -i t^/IQttI 

I inQoyj. 


output Sylvester type vectors (unmixed case) 


f indBez 


find all pure Bezout-type vectors. 


Mctke Complex 


compute tlie complex of an Tfi — vector 


ir\ "V^ T ^ n 1^ 

p±±ILT,u±OC,K.S 


pimii conipiex as tba,-'*-! a — ^ tb^-iVQ ^ 


printCohs 


prmx complex as tp-n ^u- j — r tb-n \^ ) 


IllLLJ. LIllcLU 


t^Ullo Ui Lie Vi liiajLi IJv IVl \J 1 J . OliXI r w I t/ J 


Sylvmat 


construct Sylv. matrix i^i,p — > /Co,p-i 


Bezoutmat 


construct Bezout matrix Ki a — ?• -f^o,b 


makeMatrix 


construct matrix Ki — )• 



Table 1 

The main routines of our software. 
5 Implementation 

We have implemented the search for formulae and construction of the corresponding 
resultant matrices in Maple. Our code is based on that of [7, Sect. 8] and extends 
it to the scaled case. We also introduce new features, including construction of 
the matrices of Section 4; hence we deliver a full package for multihomogeneous 
resultants, publicly available at Ww-sop. inria.fr/galaad/amantza f /soft .html [ 

Our implementation has three main parts; given data (Z, cZ, s) it discovers all possible 
determinantal formula; this part had been implemented for the unmixed case in [7]. 
Moreover, for a specific m— vector the corresponding resultant complex is computed 
and saved in memory in an efficient representation. As a final step the results of 
Section 4 are being used to output the resultant matrix coming from this complex. 
The main routines of our software are illustrated in Table 1. 

The computation of all the m— vectors can be done by searching the box defined 
in Theorem 3.4 and using the filter in Lemma 3.5. For every candidate, we check 
whether the terms K2 and K_i vanish to decide if it is determinantal. 

For a vector m, the resultant complex can be computed in an efficient data struc- 
ture that captures its combinatorial information and allows us to compute the 
corresponding matrix. More specifically, a nonzero cohomology summand K^^p is 
represented as a fist of pairs (cg,ep) where Cq = {ki,...,kt} C [l,r] such that 
q = hi = P — ^ and Cp C [0, n] with #ep = p denotes a collection of polynomi- 
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als (or a basis element in the exterior algebra). Furthermore, a term K^, is a list of 
KjjpS and a complex a list of terms K,j. 

The construction takes place block by block. We iterate over all morphisms 5a,b and 
after identifying each of them the corresponding routine constructs a Sylvester or 
Bezout block. Note that these morphisms are not contained in the representation of 
the complex, since they can be retrieved from the terms Ki^a and i^o,6- 

Example 5.1. We show how our results apply to a concrete example and demon- 
strate the use of the Maple package on it. The system we consider admits a standard 
Bezout-Dixon construction of size 6x6. But its determinant is identically zero, due 
to the sparsity of the supports, hence it neither expresses the multihomogeneous 
resultant, nor provides any information on the roots. Instead our method constructs 
a non-singular 4x4 hybrid matrix. 

Let Z = d = (1, 1) and s = (1, 1, 2). 

> read mhomo-scaled.mpl : 

> l:=vector([l,l]) : d:=l: s:= vector ( [1 , 1 ,2] ) : 

> f:= Makesystemd ,d, s) ; 

/o = «o + o,iXi + 02X2 + a^XiX2 
fi =bo + biXi + b2X2 + bsXiX2 

/2 = Co + CiXi + C2X2 + C3X1X2 + CiXi' + 05X1^X2 + 
+ CqX2 C-jX\X2 C^X\X2 

We check that this data is determinantal, using Theorem 3.7: 

> has_deter( 1, d, s) ; 

true 

Below we apply a search for all possible determinantal vectors, by examining all 
vectors in the boxes of Corollary 3.8. The condition used here is that the dimension 
of K2 and K^x is zero, which is both necessary and sufficient. 

> allDetVecsC 1, d, s) ; 

[[2, 0, 4], [0, 2, 4], [3, 0, 6], [2, 1, 6], [2, -1, 6], [1, 2, 6], [1, 1, 6], 
[1, 0, 6], [0, 3, 6], [0, 1, 6], [-1, 2, 6], [3, 1, 8], [1, 3, 8], [1, -1, 8], 
[-1,1, 8], [3, -1,10], [-1,3, 10]] 

The vectors are listed with matrix dimension as third coordinate. The search re- 
turned 17 vectors; the fact that the number of vectors is odd reveals that there 
exists a self-dual vector. The critical degree \?, p — (2, 2), thus m = (1, 1) yields the 
self-dual formula. Since the remaining 16 vectors come in dual pairs, we only men- 
tion one formula for each pair; finally, the first 3 formulae listed have a symmetric 
formula, due to the symmetries present to our data, so it suffices to list 6 distinct 
formulae. 
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Using Theorem 3.7 we can compute directly determinantal boxes: 



> detboxesC 1, d, s) ; 

[[-1,1], [1,3]], [[1,3], [-1,1]] 



Note that the determinantal vectors are exactly the vectors in these boxes. These 
intersect at m = (1, 1) which yields the self-dual formula. In this example minimum 
dimension formulae correspond to the centers of the intervals, at m = (2, 0) and 
m — (0, 2) as noted in Conj. 3.9. 



A pure Sylvester matrix comes from the vector 



> in:= vector([d[l]*convert(op(s) , ' + ')-!, -1]); 

m = (3,-1) 



We compute the complex: 



> K:= makeComplex(l,d,s,m) : 

> printBlocks(K) ; printCohs (K) ; 

Ki,2 Ko,i 

H\l, -3) © H\0, -4)2 ^ H\2, -2)2 © H\l, -3) 



The dual vector (—1,3) yields the same matrix transposed. The block type of the 
matrix is deduced by the first command, whereas printCohs returns the full de- 
scription of the complex. The dimension is given by the multihomogeneous Bezout 
bound, see Lemma 2.2, which is equal to: 



> mbezoutC 1, d, s) ; 

10 



It corresponds to a "twisted" Sylvester matrix: 



> makematrix(l,d,s,m) ; 
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-bi 


-bs 





ai 


03 














-bo 


-b2 





ao 


02 

















-bi 


-63 





ai 


as 














-bo 


-62 





ao 


ai 











-C4 


-C5 


-C8 











ai 





a3 


-ci 


-C3 


-C7 











ao 


ai 


a2 03 


-Co 


-C2 


-C6 














ao 


a2 











-C4 


-C5 


-C8 


61 





63 











-ci 


-C3 


-C7 60 bi 


62 63 











-Co 


-C2 


-cq 





60 


62 



The rest of the matrices are presented in block format; the same notation is used for 
both the map and its matrix. The dimension of these maps depend on m, which we 
omit to write. Also, B{xk) stands for the partial Bezoutian with respect to variables 

Xk. 

For m = (3, 1) we get Ki^i ® Xi,2 Ko,o, or 

H%2, 0)2 © H\0, -2)2 ^ H%3, 1) 

M(/o) 
M(/i) 

B{X2) 

Symmetric is m— (1, 3). 

For m = (3, 0), Ki^2 ^ i^o,o © KqX- 

H\l, -2) e ii"^(0, -3)2 ^ i/0(3, 0) e H\l, -2)2 








M(/o) 


B{X2) 







Symmetric is m = (0, 3). 

For m = (2, 1), we compute Ki^i ® K^^^, -f^o,o, or 

H^{1, 0)2 ® i/2(_2, -3) ^ //0(2, 1) 
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M(/i) 
_ M(/2) 

Symmetric is m = (1, 2). 
If m = (1, 1), e Xi,3 Kq^o ® Ko,2, yielding 

H%0, 0)' ® H\-3, -3) ^ 1) © H\-2, -2f 

' " 

/i 

_A(i,i),(i,i) M(/o) -M(/0_ 

We write here /j instead of M{fi), since this matrix is just the 1x4 vector of 
coefficients of /j. 

For m = (2, 0), we get Xi,2 ® i^i,3 ^ i^o,o ® ^o,i, or 

ii"^(0, -2) © H^{-2, -4) ^ i/°(2, 0) © H\0, -2) 

S(X2) 

A(2,o),(o,2) 

This is the minimum dimension determinantal complex, yielding a 4 x 4 matrix. □ 
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